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P~ ! ABSTRACT 

■ We present an analysis of three strong, intervening Mg II absorption systems (z a bs = 
| 0.603,0.656, 1.107) towards the optical afterglow of gamma-ray burst (GRB) 060418. 

From high resolution UVES spectra we measure metal column densities and find that 
' the highest redshift absorber exhibits a large amount of dust depletion compared with 

DLAs seen in QSO spectra. The intervening z a b s = 1.107 absorber is also unusual 
q ■ in exhibiting a clear 2175 A bump, the first time this feature has been definitively 

$H ' detected in a GRB spectrum. The GRB afterglow spectrum is best fit with a two 

, component extinction curve: an SMC extinction law at z = 1.49 (the redshift of the 

d • host) with E(B — V) = 0.07±0.01 and a Galactic extinction curve at z ~ 1.1 with 

E(B - V) = 0.08±0.01. We also present a moderately deep NTT i?-band image of the 
• i— i , GRB060418 field and spectroscopy of four galaxies within 1 arcminute. None of these 

1 objects has a redshift that matches any of the intervening absorbers, and we conclude 

J-j | that the galaxies responsible for the two intervening Mg II absorbers at z ~ 0.6 have 

luminosities < 0.3L*. 

Key words: quasars: absorption lines - dust, extinction - galaxies: ISM - gamma- 
rays: bursts 



1 INTRODUCTION 

One of the most powerful tools in the study of the interstel- 
lar medium (ISM) is the use of bright background sources 
on whose continua the absorption signature of intervening 
gas is clearly scribed. In the Milky Way, these background 
sources are usually stars and for extragalactic work, quasars. 
Combined with the large apertures of modern telescopes and 
the high resolving power of echelle spectrographs, absorp- 
tion line spectroscopy has yielded a detailed model of the 
ISM and the intergalactic medium (IGM) from the present 
day to z—5 (see reviews by Savage & Sembach 1996; Rauch 
1998; Wolfe, Gawiser & Prochaska 2005). Very recently, a 
new type of background source has become available - the 



optical afterglows of gamma-ray bursts (GRBs). In partic- 
ular, the launch of the Swift satellite and its rapid provi- 
sion of accurate coordinates, has seen the coming of age of 
prompt, optical, ground-based follow-up of GRBs. Although 
low resolution spectra of the optical afterglow obtained a few 
hours after the GRB already showed signs of intervening (as 
well as host galaxy) absorption (e.g. Metzger et al. 1997; 
Mirabal et al. 2002; Vreeswijk et al. 2004), a water-shed 
has been recently crossed with the first echelle spectroscopy. 
In a number of cases, this follow-up is obtained within an 
hour of the burst (e.g. Chen et al. 2005), with the fastest 
follow-up being a mere 10 minutes after the GRB trigger 
(Vreeswijk et al., in preparation). A number of exciting in- 
sights into the GRB environment are emerging from these 
prompt, high resolution data, including the detection of nu- 
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Figure 1. i?-band image of the field of GRB060418 taken with 
EMMI on the NTT. The objects for which we obtained spectra are 
marked (A-D), as is the approximate position of the GRB (cross). 
Galaxies A-D have the following i?-band magnitudes and red- 
shifts (based on absorption line matches with an elliptical galaxy 
template): A - R = 21.6, z = 0.22; B - R = 21.4, z = 0.29; C 
- R = 21.1, z = 0.36; D - _R = 21.1, z = 0.34. Galaxy E has a 
magnitude of R = 22.7. 



2 DATA 

2.1 UVES spectroscopy of the afterglow 

Our group has been involved in commisioning a Rapid Re- 
sponse Mode (RRM) 1 at the Very Large Telescope (VLT) in 
Chile. In this paper, we discuss the intervening absorption 
systems of GRB 060418 which was observed with the RRM 
using the UVES echelle spectrograph. Details of the obser- 
vations can be found in Vreeswijk et al.(in preparation). In 
brief, the wavelength coverage is 330 - 1000 nm, the reso- 
lution ~ 7 km s _1 and the S/N per pixel ranges from 15 - 
25. Special care was taken with the merging of the echelle 
orders, combining the coverage from UVES's three CCDs 
and flux calibration (which included a correction for Galac- 
tic extinction with E(B — V) = 0.224, Schlegel, Finkbeiner 
& Davis 1998). We selected strong Mg II absorbers with rest 
frame EW(Mg II A 2796) > 1 A; these are listed in Table 
In none of these Mg II absorbers do we cover the Lya 
line, but we can apply the Ellison (2006) D-index as an in- 
dicator of whether the absorber is likely to be a DLA. In 
their pilot sample, Ellison (2006) found that absorbers with 
D = 1000 x EW(Mg II A 2796) / Av > 6.3 have a 90% 
probability of being DLAs. All three Mg II systems in our 
sample have U-indices greater than 6.3. 



merous fine structure transition lines (Vreeswijk et al. 2004; 
Chen et al. 2005; Prochaska, Chen & Bloom 2006). 

With a growing database of afterglow spectra, and with 
the addition of high S/N, high resolution spectra in a num- 
ber of cases, we are also approaching the point at which 
preliminary statistics of intervening absorbers can be at- 
tempted. Prochter et al. (2006) have recently reported a puz- 
zling over-abundance of intervening high equivalent width 
(EW) Mg II absorbers compared with QSO sightlines. Three 
possible explanations are proposed: dust extinction in QSO 
sightlines (see also Savaglio, Fall & Fiore 2003), gravitational 
lensing in GRB sightlines and associated (intrinsic) absorp- 
tion as responsible for some of the GRB Mg II population. 
Although the first of these, bias due to dust, has been shown 
to be small towards QSOs at absorption redshifts spanning 
0.6 < z abs < 3.5 (Ellison et al. 2001, 2004; Akerman et 
al. 2004; Ellison, Hall & Lira 2005; Jorgenson et al. 2006), 
some concerns about sample size remain and GRBs offer an 
interesting alternative probe. 

Another attractive property of studying intervening ab- 
sorbers towards GRB afterglows is that the fading source al- 
lows for a more sensitive search for absorber-galaxy counter- 
parts than is possible towards QSOs. Although the canonical 
result for Mg II absorbers is that they are associated with 
relatively bright galaxies (e.g. Steidel, Dickinson & Persson 
1994) there are issues concerning incompleteness and mis- 
identification (Churchill, Kacprzak & Steidel 2005). In this 
paper, we combine high resolution optical spectroscopy of 
the afterglow to determine dust and abundance properties 
of the intervening absorbers, with imaging of the GRB field 
and spectroscopy of low impact parameter galaxies. 

We use the following cosmological parameters: Qm = 
0.3, ft A = 0.7, Ho =70 km/s/Mpc. 



2.2 Imaging 

We obtained 2 x 300 second exposures of the GRB field on 21 
April 2006 using EMMI-RILD on the NTT at La Silla. The 
final image is shown in FigureQand has a 5a detection limit 
of R = 23.45 within a 3 arcsecond diameter aperture. We 
also obtained spectra of four objects around the host galaxy 
(galaxies A, B, C and D), with impact parameters ranging 
from 13 to 60 arcseconds 3 . The spectra were obtained with 
the EMMI grism #3 yielding a FWHM resolution of 8 A. 
The total integration time per object was 7200 seconds for 
galaxies A and B and 2400 seconds for galaxies C and D. 
None of the 4 objects exhibit strong emission lines but a 
cross-correlation with an elliptical galaxy template yields 
redshifts below 0.4 (see Figure caption) based on stellar 
absorption features. None of these objects is therefore likely 
to be associated with the strong absorption systems listed 
in Table The second closest object (after galaxy A) to 
the GRB that is detected in our image is galaxy E with an 
impact parameter of 23 arcseconds and an approximate R 
band magnitude of 22.7, too faint to follow-up with a 4-m 
class telescope. 
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Table 1. Metal Column Densities of Intervening Absorption Systems 





EW(Mg II A 2796) rest 


D-index 


N(Fe II) 


N(Si II) 


N(Zn II) 


N(Cr II) 


N(Mn II) 


N(Ti II) 


0.603 


1.24±0.09 A 


8.0 


14.85±0.20 








12.79±0.03 


12.51±0.03 


0.656 


1.00±0.08 A 


7.2 


13.99±0.03 








<12.15 




1.107 


1.87±0.06 A 


7.1 


14.67±0.03 


15.47±0.05 


12.87±0.03 


< 12.70 


13.04±0.06 


< 11.75 



All limits are 5o\ 
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Figure 2. Example Voigt profile fits for the z^bs = 1.107 ab- 
sorber towards GRB060418. 



derlying nucleosynthetic patterns plus the relative suscepti- 
bility to dust depletion. Disentangling these effects in order 
to profile the chemical enrichment (and, therefore, star for- 
mation) histories has been a long-term goal of DLA research 
(e.g Dessauges-Zavadsky et al. 2006). We do not know the 
absolute metallicity of the z a b s = 1.107 GRB060418 absorber 
due to the lack of N(HI) measurement. However, the large 
depletion factor (see below) implies that the metallicity may 
be at least 1/10 of solar (Ledoux et al. 2003). Indeed, this 
is the average metallicity of DLAs at z ~ 1 (Meiring et al. 
2006), which would require a relatively high N(HI)=21.24. 
Alternatively (and taking Zn as an undepleted measure of 
metallicity), this absorber would have solar metallicity for 
N(HI) = 20.24. 

Even in the absence of absolute abundances, the relative 
abundances that can be deduced from Table reveal that 
the 2 a bs = 1.107 GRB060418 absorber has a large amount 
of dust depletion (i.e. low gas-to-dust ratio). This is usually 
inferred from the ratios of iron or chromium to zinc (e.g. Pet- 
tini et al. 1997), relative to the solar value: [Fe/Zn], [Cr/Zn]. 
We determine [Fe/Zn] = — 1.04 which is lower than ~ 95% of 
DLAs (e.g. Meiring et al. 2006) 5 . We also measure [Cr/Zn] 
< —1.19 which is lower than all of the DLAs in the compila- 
tion of ~ 40 DLAs presented by Meiring et al. (2006). Other 
abundance ratios show similar extreme values, also indica- 
tive of dust. For example Si and Ti are both alpha capture 
elements whose Galactic stellar abundances track each quite 
well, although with a ~ 0.1 dex overabundance of Si (e.g. 
Fulbright 2000). However, the abundances of Ti and Si are 
very discrepant here: [Ti/Si] < —1.10. Again, this is likely 
due to the highly refractory nature of Ti. 



3 RELATIVE ABUNDANCES IN GRB 
ABSORBERS 

We have fitted the metal line profiles of intervening ab- 
sorbers in the UVES spectrum with Voigt profiles using the 
software VPFIT . The column densities are listed in Ta- 
blc0 Unfortunately, the two lower redshift systems (z a bs = 
0.656, 0.603) do not have coverage of many metal lines, al- 
though the higher redshift z a b s = 1.107 has numerous tran- 
sitions included in our wavelength coverage; a selection of 
the absorption lines and fits is shown in Figure 

Abundances will be driven by the combination of un- 



1 http:/ /www. eso.org/observing/p2pp/rrm. html 

2 Note that the EWs quoted for Mg II are slightly higher than 
the values used for the .D-index calculation. In Table HI we quote 
the full line EW, but the D-index excludes the extreme velocities 
of the profile where the absorption has a significance of < 3cr 

3 1 arcsecond is 6.69 /i~~g kpc (proper) at z = 0.6 and 8.17 h^ 
kpc at 2 = 1.1. 

4 http:/ /www. ast.cam.ac.uk/~rfc/vpfit. html 



4 EXTINCTION CURVE 

In order to quantify the amount of extinction that is asso- 
ciated with the ;z a bs = 1.107 absorber, we investigate fitting 
the flux calibrated spectrum as a power law continuum com- 
bined with extinction curves for the Milky Way, SMC and 
LMC as parameterised by Pei (1992). The power-law index 
(/9) and E(B — V) are allowed to vary as free parameters. 
We included two dust components, one fixed at the redshift 
of the host (z — 1.49) and one for an intervening absorber 
whose redshift is also a free parameter. In the absence of a 
2175 A feature, the GRB host galaxy extinction is best fit 
with an SMC curve. However, at z ~ 1.1 we find a clear 2175 
A bump, indicative of carbonaceous grains, which requires 
either an LMC or MW extinction curve to achieve a good 
fit. We adopt the Galactic extinction curve as the best fit, 

5 Where necessary, we adopt the solar abundance scale of Lodders 
(2003). 
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5 DISCUSSION AND CONCLUSIONS 

We have presented the abundances of intervening absorbers 
towards GRB060418 and, for the first time in a GRB after- 
glow spectrum, detected the 2175 A bump associated with 
an absorber at z a b s = 1.107. The GRB afterglow spectrum 
can be well fitted with a Milky Way extinction curve with 
Zdust ~ 1.1 and E(B - V) = 0.08, combined with SMC dust 
at the redshift of the host. The presence of dust in this ab- 
sorber is supported by significantly sub-solar [Fe/Zn] and 
[Cr/Zn] ratios, which are similar to the highest dust-to-gas 
ratios measured in QSO DLAs. 

It has been previously suggested (e.g. Savaglio et al. 
2003) that GRBs may reveal dustier absorbers than QSO 
surveys, since bright afterglows could be observed even in 
the presence of substantial foreground extinction. Recently, 
Vladilo & Peroux (2005) showed that there is a relation- 
ship between N(Zn) and Ay in the Galaxy and in Vladilo et 
al. (2006) suggested a widespread correspondence between 
depletion and reddening. The N(Zn) - Ay relation applied 
to the z a bs = 1.107 absorber towards GRB060418 predicts 
a rest-frame Ay = 0.21. This is in good agreement with 
the value determined from the Galactic extinction curve fit 
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based on a lower x 2 , which we show in Figure [3] The best 
fit parameters for intervening Galactic-type dust are (rest- 
frame, with la errors quoted) E(B - V) = 0.08±0.01 {A v 
= 0.25) at z = l.llSiHoi and for the SMC host extinc- 
tion, E(B - V) = 0.07±0.01 (Av = 0.22). The underlying ^ 1.5 - 
power- law continuum has a best-fit value of (3 = — 0.291q'q*, 
again with la errors quoted based on 68% confidence limits 
from the x 2 minimization. The relatively low reddening at 
z ~ 1.1 is consistent with the weak Ca II absorption, which 
we measure to have a rest frame equivalent width of only 0.1 
A (Wild, Hewett & Pettini 2006). The extinction in the host 
galaxy is consistent with the (rather wide) range of values 
observed in other GRB afterglow spectra (e.g. Galama & 
Wijers 2001; Vreeswijk et al. 2004). 

Given the large Galactic reddening in the direction of 
the GRB, we repeat the fitting process with different values 
of z — extinction in order to assess the sensitivity of the 
fit to this parameter. Changing the Galactic E(B — V) by 
0.1 magnitudes results in a negligible effect on the best fit 
E(B - V) of the z ~ 1.1 absorber (A E(B - V) = 0.003), 
although the best fit /3 changes by up to A/3 = 0.2. 

Our data provide the first definitive detection of the 
21 75 A bump in a GRB spectrum and one of the highest red- 
shift detections of the feature in the literature. The bump has 
only previously been observed in a handful of QSO absorbers 
(e.g. Wang et al. 2004; Junkarinnen et al. 2004) , some gravi- 
tational lenses (e.g. Motta et al. 2002; Wucknitz et al. 2003) 
plus a tentative detection associated with the host galaxy 
of GRB991216 (Vreeswijk et al. 2006). Reddening in QSO 
absorbers is usually best fit by an SMC-type extinction (e.g. 
York D. G. et al. 2006). One interesting feature of our de- 
tection is that if the bump occurs at z = 1.107 (the redshift 
of the observed metal lines), its inferred rest wavelength is 
2186 A. This is within the ± 17 A scatter seen in Galactic 
sightlines, despite the theoretical prediction that more vari- 
ation should be seen for a graphite carrier (e.g. Fitzpatrick d 

& Massa 1986). & 

a 
-1 
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Figure 3. Flux calibrated spectra and errors (black lines) and 
fits (red) show the intrinsic power law spectrum (upper curve) 
and dust reddened spectrum (lower curve). The raw spectrum 
has been smoothed with a Gaussian of FWHM=0.25 A for pre- 
sentation purposes. A two-component dust fit is used with an 
SMC extinction curve fixed at z = 1.49 for the host and an inter- 
vening source of dust (redshift is fitted as a free parameter) with 
a MW extinction curve. 
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Figure 4. The observed frame S-band extinction induced by 
3 different dust models as a function of the absorber (dust) and 
redshift. The dot-dash curve represents the dust present in the 
•2abs = 0.524 DLA towards the QSO AO 0235+164, the dustiest 
DLA known. The dashed curve illustrates extinction for the best 
fits for the z a b s ~ 1.1 absorber seen towards GRB060418. The 
solid curve represents the maximum extinction inferred from the 
optically complete DLA survey of Ellison, Hall &: Lira (2005). 



(A v = 0.25). As pointed out by Vladilo & Peroux (2005), 
this amount of extinction is quite modest, and unlikely to 
produce a significant bias in absorber statistics. We demon- 
strate this in Figure ^] where we show the observed-frame 
extinction in the B-band of a z = 3 QSO as a function of 
Zdust- We calculate the change in magnitude due to extinc- 
tion by taking the composite SDSS QSO spectrum (Vanden 
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Berk et al. 2001), applying a typical B-band transmission 
filter (we use the filter curve from SUSI2 on the NTT) and 
applying the extinction curves of Pei (1992). We show three 
curves, illustrating the case for the z a b s = 1.107 absorber to- 
wards GRB060418 and also the z abs = 0.524 DLA towards 
AO 0235+164. The latter is one of the dustiest DLAs known 
with E(_B — V) = 0.23, close to solar metallicity and even 
exhibiting diffuse interstellar bands (Junkarinnen et al.2004; 
York, B. A. et al. 2006). In addition, we show the effect of 
SMC extinction when E(B - V) = 0.04, the upper limit de- 
termined by Ellison et al. (2005) from B — K colours of an 
optically complete QSO sample. Only the model based on 
the DLA towards AO 0235+164 introduces more than ~ 0.7 
mags of extinction in the observed frame B-band. Moreover, 
although the z a b s = 1.107 absorber towards GRB060418 has 
extreme abundance ratios and exhibits a 2175 A bump, the 
observed extinction is comparable to that determined from 
complete optical QSO surveys (i.e. the results of Ellison et 
al. 2005). Even with dust properties as unusual as this (i.e. 
extreme depletions and a 2175 A bump), we support the 
conclusion of Prochter et al. (2006) that there is currently 
no convincing evidence to suggest that dust bias is the cause 
of an enhanced Mg II number density (see also Frank et al. 
2006). 

Our image of the field around GRB 060418 reveals an 
absence of bright galaxies (R < 23.45) close to the host 
galaxy. The faintest galaxy within 30 arcsecs of the after- 
glow (galaxy E on Figure would have a luminosity of 0.7 
L* at z = 0.65 6 . Although this is consistent with the lumi- 
nosities of Mg II absorbers found by Steidel et al. (1994), 
the impact parameter is ~ 200 h^g kpc, so we conclude 
that this is unlikely to cause strong Mg II absorption in the 
afterglow spectrum. Our limiting magnitude of R > 23.45 
corresponds to ~ 0.3L* at z — 0.65 indicating that the ab- 
sorbing galaxies at z ~ 0.6 are, in this case, fainter than 
many of the galaxies which have been associated with QSO 
Mg II absorbers. For example, even when the glare of the 
quasar was not a contaminating factor, O'Meara, Chen & 
Kaplan (2006) found ~ L* galaxies associated with QSO 
Mg II absorbers. 

In summary, based on this one field, we have found one 
absorber that is apparently unusually depleted compared 
with QSO DLAs and is rare in exhibiting a Galactic extinc- 
tion curve with a 2175 A bump. This is the first definitive 
detection of the 2175 A dust bump due to a galaxy in the line 
of sight to a GRB afterglow. We have also shown that the 
galaxies responsible for two further absorbers at z ~ 0.6 are 
fainter than 0.3L*. Determining the extinction properties 
of more intervening systems in GRB sightlines is required 
to determine how common these properties are in galaxies 
seen in absorption in GRB afterglow spectra. The combina- 
tion of prompt and positionally accurate Swift triggers with 
ground-based follow-up by teams such as our group 7 and 
GRAASP 8 is likely to provide significant progress in this 
direction in the near future. 

6 We used the value of L* at 0.6 < z < 0.8 from Ilbert et al. 
(2005), including a correction to the Vega system and assuming 
that the observed frame R is similar to the rest frame B band at 
z ~ 0.6. 

7 http:/ /www. sc. eso.org/~grbalert/uves/home. html 

8 http:/ /graasp. ucolick.org/graasp. html 
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